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Abstract 

We investigate deviation from the tri-bimaximal mixing in the case of ultra 
high energy neutrino using ICECUBE detector. We consider the ratio of 
number of muon tracks to the shower generated due to eletrons and hadrons. 
Our analysis shows that for tri-bimaximal mixing the ratio comes out around 
4.05. Keeping 612 and 623 fixed at tri-bimaximal value, we have varied the 
angle 6^13 = 3°, 6°, 9° and the value of the ratio gradually decreases. The vari- 
ation of ratio lies within 8% to 18% from the tri-bimaximal mixing value and 
it is very difficult to detect such small variation by the ICECUBE detector. 
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1 Introduction 



Various experiments for solar and atmospheric neutrinos provide a range for 
the values of solar mixing angle ^0 = 612 (the 1-2 mixing angle) [1] that 
corresponds to solar neutrino oscillations and also a range for atmospheric 
mixing angle ^atm = ^23 (the 2-3 mixing angle) [1] around their best fit values. 
The tri-bimaximal mixing condition of neutrinos are given by sin ^12 = 
sin ^23 = and sin ^13 = [2]. Possible deviations from tri-bimaximal 
mixing can be obtained by probing the ranges of 9i2 and ^23 given by the 
experiments. Also the exact 13 mixing angle ^13 is not known except that 
the CHOOZ [3] gives an upper limit for ^13 (< 9°). Probing the deviations of 
012 and ^23 for different values for ^13 is significant not only to understand 
the neutrino flavour oscillations in general but also for the purpose of model 
building for neutrino mass matrices. 

In this work wc explore the possibility for ultra high energy (UHE) neu- 
trinos from distant Gamma Ray Bursts (GRBs) for probing the signatures 
of these deviations of the values of the mixing angles from tri-bimaximal 
mixing as discussed above. One such proposition of using UHE neutrinos 
is described in a recent work by Xing [4]. Gamma Ray Bursts are short 
lived but intense burst of gamma rays. During its occurence it outshines all 
other luminous objects in the sky. Although the exact machanism of GRBs 
could not be ascertained so far but the general wisdom is that it is pow- 
ered by a central engine provided by a failed star or supernova that possibly 
turned into a black hole, accretes mass at its surroundings. This infalling 
mass due to gravity bounces back from the surface of black hole much the 
same way as the supernova explosion mechanism and a shock is generated 
that flows radially outwards with enormous amount of energies (~ 10^^ ergs). 
This highly energetic shock wave drives the mass outwards, in the form of 
a "fireball" that carries in it, protons, 7 etc. The pions are produced when 
the accelerated protons inside the fireball interacts with 7 through a cosmic 
beam dump process. UHE neutrinos are produced by the decay of these 
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pions. Thus a generic cosmic accelerator accelerates the protons into very 
high energies which then beam dump on 7 in the "fireball" as also at the 
cosmic microwave background (CMB) and ultra high energy neutrinos are 
produced. 

The GRB neutrinos, due to their origin at astronomical distances from 
earth, provide a very long baseline for the earth bound detectors for UHE 
neutrinos such as ICECUBE [5]. The oscillatory part of the neutrino flavour 
oscillation probabilities (sin^(Am^[L/4E])) averages out to 1/2 becuase of 
this very long baseline L (~ hundreds of Mpc) and the Am^ (mass square 
difference of two neutrinos) range obtained from solar and atmospheric neu- 
trino experiments are Am2i ~ 10~^ eV^ and Am32 ~ 10"^ eV^ respectively 
(L/Am^ >> 1). Thus for neutrino flavour oscillation, in this case, the effect 
of Am^ is washed out and governed only by the three mixing angles namely 
^12 = ^0, ^23 = ^atm ^nd ^13. The purpose of the present work is to probe 
whether or not the possible variations of 612 and ^23 from their best fit values 
can be ascertained by UHE from distant GRBs. 

The GRB neutrinos, on arriving the earth, undergo charged current (CC) 
and neutral current (NC) interactions with the earth rock and the detec- 
tor material. The CC interactions of produce secondary muons and the 
same for electrons produce electromagnetic shower (i/^ -\- N ^ /i + X and 
I'e + N ^ e + X) . The former will produce secondary muon tracks and can 
be detected by track-signal produced by the Cerenkov light emitted by these 
muons during their passage through a large underground water/ice Cerenkov 
detectors like ICECUBE. The ICECUBE is a Ikm^ detector in south pole 
ice and can be considered to be immersed in the target material for the UHE 
neutrinos where the neutrino interactions are initiated. In case of z/g, the 
electrons from the UgN CC interactions, shower quickly and can also be de- 
tected by such ICECUBE detector. The case of Ut is somewhat complicated. 
The first CC interaction of Ur {vt + N ^ t + X) produces a shower ("first 
bang") alongwith a r track. But the Vr is regenerated (with diminished en- 
ergy) by the decay of r and in the process produces another hadronic or 
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electromagnetic shower ("second bang"). The whole process is called double 
bang event. In case the first bang could not be detected, then by possible 
detection of second bang (with showers) the r track can be reconstructed or 
identified and this scenario (the r track and the second bang) is called the 
lollipop events. An inverted lollipop event is one where only the first bang 
{ut- + N —>■ T + X) is detected and the subsequent r track is detected or 
reconstructed. As mentioned in Ref. [6], the detection of from their CC 
interaction mentioned above is not every efficient by a Ikm^ detector since 
the double bang events can possibly be detected only for the Ut energies be- 
tween 1 PeV to 20 PeV beyond which the tau decay length is longer than 
the width of such detector and at still higher energies the flux is too small 
for such detectors for their detection. Hence, in the present work we do not 
consider the events initiated by Ut-N CC interactions. However, for i/j- we 
consider the process that may yield events higher than the "double bang" 
events. We consider the decay channel of r lepton [7], obtained from charged 
current interactions of i/r, where muons are produced (i/,- — > r — > Dfj^^Vr) 
which can then be detected as muon tracks [8] in ICECUBE detector. The 
neutral current (NC) interactions of all flavours however will produce the 
shower events at ICECUBE and they are considered in this investigation. 

This paper is organised as follows. In Section 2 we describe the formahsm 
for neutrino fluxes of the three species while reaching the earth. The nature 
of the GRB flux taken for present calculations is also discussed. The flux 
suffers flavour oscillations while traversing from GRB site to the earth. The 
oscillation probabilities are also calculated and the oscillated flux obtained 
on reaching the earth is determined. They are given in Section 2.1. We also 
describe in this section the analytical expressions for the yield of secondary 
muons and shower events at the ice Cerenkov kilometre square detector like 
ICECUBE. This is given in Section 2.2. The actual calculations and re- 
sults are discussed in Section 3. Finally, in Section 4, some discussions and 
summary are given. 
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2 Formalism 



2.1 GRB Neutrinos Fluxes 

The neutrino production in GRB is initiated through the process of cosmo- 
logical beam dump by which a highly accelerated protons from GRB interacts 
with 7 to produce pions which in turn decays to produce ^'^(tv) ^e(i^e) 
much the same ways as atmospheric neutrinos are produced. They are pro- 
duced in the proportion 2v^ : 2z7^ : lu^ : li7e [9]. 

For the present calculation we consider the isotropic flux [10] resulting 
from the summation over the sources and as given in Gandhi et al [11]. The 
isotropic GRB flux for i/^ + is given as 



In the above, 

= 4.0 X 10-^^ n=l, for < 10^ GeV 
7\A = 4.0 X 10-^ n = 2, for > 10^ GeV 



Thus, 



dN dN- 



(2) 



dE^ dE^ 
The neutrinos undergo flavour oscillation during their passage from the GRB 
to the earth. Under three flavour oscillation, the v^. and originally created 
at GRB will be oscillated to Vr. Thus after flavour oscillations, the v^. fluxes 
(F^J, i/^ fluxes (-F^^), fluxes (F^^) become 

Fur = Pue^UrfpUe + Puu.^Ur'i^U^ ■ (3) 
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The transition probability of a neutrino of flavour a to a flavour (3 is given 



by, 



5ap -^^ U^^UpV^Vp^ sin' (^^j (4) 



In the above oscillation length Aj, is given by 



(5) 



Because of astronomical basehne Am'^L/E » 1, the oscillatory part be- 
comes averaged to half. Thus, 



Therefore 



sm 



(6) 



<^a/3 - 2 X! UaiUp.UajUp. 



j>i 



i 



(7) 



where use has been made of the condition Ylii VaJJ^i = ^a^- 
With Eq. (7), Eq. (3) can be rewritten in matrix form 
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In Eq. (8) above, we have used the initial flux ratio from GRB to be 0, 
4>u^ : 4>ur = 1:2:0. From Eq. (8) it then follows that, 

F.. = {u'ei[i + (ul,-u^,)] + uUi + (u%-u^,)]+ 



F., = K'l[l + (f/^\-f/.\)]+f/i[l + (^i-^r2)] + 



IJ.2 



t^r3[l + (t^i-t^r3)]}0.e 

The MNS mixing matrix U for 3-flavour case is given as 



U 



C12C13 S12C13 
-C23S12 - S23S13C12 C23C12 - S23S13S12 

V S23S12 — C23S13C12 — •S23C12 ~ C23S13S12 



Sl3 \ 
S23C13 
C23C13 / 



(9) 



(10) 



We are not considering any CP violation here. Hence Eqs. (3) - (9) above 
also hold for antineutrinos. 



2.2 Detection of GRB neutrinos 

The z/^'s from a GRB can be detected from the tracks of the secondary muons 
produced through the CC interactions. 

The total number of secondary muons induced by GRB neutrinos at a 
detector of unit area is given by (following [12, 9, 13]) 

S= dE,-^P,,,,(E,)P^(E,,E,^r) (11) 

In the above, Pgurv is the probability that a neutrino reaches the detector 
without being absorbed by the earth. This is a function of the neutrino- 
nucleon interaction length in the earth and the effective path length X{6z) 
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(gm cm~^) for incident neutrino zenith angle 9^ {9z = for vertically down- 
ward entry with respect to the detector). This attenuation of neutrinos due 
to passage through the earth is referred to as shadow factor. For an isotropic 
distribution of flux, this shadow factor (for upward going neutrinos) is given 

by 

P.urv{E.) = ^ /_° d cose I d0exp[-X(^,)/Lint]. (12) 
where interaction length Lint is given by 

^'■^^ ^ a^^ik)NA ^^^^ 

In the above iV^(= 6.022 x lO^^gm"^) is the Avogadro number and o'*°*(= 
o"^*" + (7^^) is the total cross section. The effective path length X{9z) is 
calculated as 

X{e,) = J p{r{0,,i))di. (14) 

In Eq. (9), p{r{9z,i) is the matter density inside the earth at a distance 
r from the centre of the earth for neutrino path length i entering into the 
earth with a zenith angle Oz- The quantity Pf^{E^, Ethr) in Eq. (6) is the 
probability that a secondary muon is produced by CC interaction of z/^ and 
reach the detector above the threshold energy i^thr- This is then a function 
of u^N (N represents nucleon) - CC interaction cross section cr*"^ and the 
range of the muon inside the rock. 

P^{E,, Eti^r) = Na<J^^{R{E,; £;thr)) (15) 

In the above {R{Ei,; £^thr)) is the average muon range given by 

{R{E,; EthO) = ^ / dyR{E,{l - y), ^thr) ) "'^^ (16) 

(7^^ Jo ay 

where y = [E^ — E^)/Ey is the fraction of energy loss by a neutrino of energy 
El, in the charged current production of a secondary muon of energy E^. 
Needless to say that a muon thus produced from a neutrino with energy 
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El, can have the detectable energy range between E'thr and E^,. The range 
R{En, £^thr) for a muon of energy E"^ is given as 

The average lepton energy loss with energy E^ per unit distance travelled is 
given by [12] 

The values of a and (3 used in the present calculations are 

a = {2.033 + 0.077 ln[E^(GeV)] } X lO^^GcVcm^gm-^ 

(3 = {2.033 + 0.077 ln[£;^(GeV)] } X 10" Wgm"^ (19) 

for E^ ^ 10^ GeV [14] and 

a = 2.033 X lO-^GeVcm^gm-^ 

P = 3.9 X lO-^cm^gm-^ (20) 

otherwise [15]. For muon events obtained from i/^ CC interactions, in 
Eq. (11) will be replaced by F^,^ (Eq. 9). 

As discussed earlier, the events due to Ut CC interactions is considered 
only for the process where the decay of secondary r lepton produces muon 
which then detected by the muon track. The probability of production of 
muons in the decay channel r v^^Vr is 0.18 [7, 8]. The generated muon 
carries a fraction 0.3 of energy of original Vr (a fraction 0.75 of the energy of 
the is carried by secondary r lepton and a fraction of 0.4 of r lepton energy 
is carried by the muon [7, 12, 8]). For the detection of such muons, the Eqs. 
(10 - 16) is applicable with properly incorporating the muon energy described 
above. Needless to say, in this case, in Eq. (11) is to be replaced by Fj,^ 
(Eq. 9). 
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For the case of showers, we do not have the advantage of a specific track 
and then the whole detector volume is to be considered. The event rate for 
the shower case is given by 

N,^ = J dE,^P,,„{E,) X J lj^P^,{E,,y) . (21) 

In the above, — a^^ (for electromagnetic shower from i/g charged current 
interactions) or a^'^ as the case may be. In the above Pint is the probability 
that a shower produced by the neutrino interactions will be detected and is 
given by 

Pint = pNact'L (22) 

where p is the density of the detector material and L is the length of the 
detector (L = 1 Km for ICECUBE). 

For each case of shower events, in Eq. (21) is to be replaced by F,y^ 
or F,y^ or F^,^ as the case may be. 

3 Calculations and Results 

The secondary muon yield at a kilometre scale detector such as ICECUBE 
is calculated using Eqs. (6 - 20). The earth matter density in Eq. (9) 
is taken from [9] following the Preliminary Earth Reference Model (PREM). 
The interaction cross-sections - both charged current and total - used in these 
equations are taken from the tabulated values (and the analytical form) given 
in Ref. [11]. In the present calculations E^^raax — 10^^ GeV and threshold 
energy E'thr = 1 TeV are considered. 

For our investigations, we first define a ratio TZ of the muon events (both 
from i/^ (and P^^) and Uj- (and P,-)) and the shower events. As described in 
the previous sections, the muon events are from i/^ (and P^^) and Uj- (and Pr), 
whereas the shower events include electromagnetic shower initiated by CC 
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interaction of v^. and NC interactions of neutrinos of all flavours. Therefore, 

^ = ^ (23) 

where, 

= S{ioY v^) + 5(for vr) 
^sh = -^sh(for i/e CC interaction) 
+ A^sh(for fe NC interaction) 
+ A^sh(for NC interaction) 
+ A^sh(for Vr NC interaction) (24) 

The purpose of this work is to explore whether UHE neutrinos from ORB 
will be able to distinguish any variation of Q\2 and Q2?, from their best fit 
values. The tri-bimaximal mixing condition is denoted by the best fit values 
of Qyi and ^23 for ^13 = 0". The best fit value of Qyi = 35.2" and that of 

= 45°. We first vary Ou in the hmit 30° < 6^12 < 38° and vary 6*23 in the 
limit 38" < 612 < 54° with ^13 = and for each case calculate the ratio TZ 
using Eqs. (1 - 24). We find that TZ varies from 3.14 to 4.25. One readily sees 
that the variation in muon to shower ratio is not very significant. The fiux 
and other uncertainties of the detector may wash away this small variations. 
TZ obtained from tri-bimaximal condition given above is 4.05. 

The same operation is repeated for three different values of 613, namely 
612 = 3° 6° and 9° with similar results. The results are tabulated below. 



^13 


"^Max 


TZMin 


7^ at 

612 = 35.2°, ^23 = 45° 


0° 


4.78 


3.80 


4.05 


3° 


4.75 


3.77 


4.01 


6° 


4.72 


3.75 


3.98 


9° 


4.69 


3.73 


3.96 
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Table 1. Maximum and minimum values of ratio TZ for different values of mixing 
angles 

We have also plotted the variation of TZ with 612 and 623 for four fixed 
values of ^^13 as given in Table 1. These are shown in Figs la - Id for 613 = 
0°, 3° , 6° and 9° respectively. 

As is evident from Table 1 and Fig. 1, the variation of muon tracks 
to shower ratio is not very significant with the deviation from the best fit 
values of the mixing angles. The ratio TZ varies upto only ~ 18%. We 
have also calculated the muon track signal for 1 year of ICECUBE run. For 
^13 — 0, this varies from ~ 99 to ~ 115, whereas the muon yield obtained 
for tri-bimaximal mixing is 103. So the variation for deviation from tri- 
bimaximal mixing condition is between 4% — 11%. This variation is also not 
significant given the sources of uncertainty in the flux and the sensitivity of 
the ICECUBE detector. Firstly, the flux itself can be uncertain by several 
factors. This can induce errors in calculation of muon yield and shower rate. 
If the flux uncertainties are energy-dependent, even the ratio TZ can also be 
affected. Also the simulation results for ICECUBE detector by Ahrens et al 
[16] shows the cosmic neutrino signal is well below the atmospheric neutrino 
background for one year data sample after applying suitable cuts (for the 
source flux x dN^/dE^ = 10~'''cm^s~^sr~^GeV). The diffuse flux needed 
for a 5(j significance detection after 1 year is well below the experimental 
limits [16, 17]. There can also be systematic uncertainty arises out of optical 
module (OM) sensitivity which is affected by the rcfrozcn ice around OM, 
optical properties of the surrounding ice, trapped air bubbles in the OM 
neighbourhood etc. An estimation of these uncertainties for a signal is 
calculated to be around 20% [16]. Taking into account these uncertainties 
and sensitivity limit, it is difficult by a detector like ICECUBE to detect the 
deviation 18%), if any, from tri-bimaximal mixing through the detection 
of UHE neutrinos from a GRB. 
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4 Summary and Discussions 



In summary, we investigate the deviation from the well known tri-bimaximal 
mixing in the case of Ultra High Energy neutrinos from a Gamma Ray Burst 
detected in a kilometer scale detector such as ICECUBE. We have calculated 
the ratio TZ of the muon track events and shower events (electromagnetic 
shower from charged current interactions of ly^ and hadronic showers from 
neutral current interactions of neutrinos of all flavours) for tri-bimaximal 
mixing condition given by 9^ = 35.2°, 623 = 45.0°, ^13 = 0°. We then 
investigate the possible variation of TZ from tri-bimaximal mixing condition 
by varying $12 and $23 within their experimentally obtained range for four 
different values of ^13 namely 0°, 3° , 6° and 9°. 

The isotropic flux of GRB neutrinos are obtained following Waxman- 
Bahcall [10] type parametrization of the flux and summation over the sources. 
The initial parametrization of neutrino flux can be written as 



where is the spectral break energy (~ 10^ GeV) and is related to photon 
spectral break energy, Lorentz factor etc. 

The GRB neutrinos after reaching the earth has to pass through the earth 
rock (for upward going events) to reach the detector to produce muon tracks 
or shower. In the calculation therefore, the attenuation of neutrinos through 
the earth (shadow factor) is estimated. The muons produced out of charged 
current interactions of neutrinos should also survive to enter the detector 
and produce tracks. Therefore, to estimate the muon track events, the en- 
ergy loss of muons through the rock is also estimated. The average lepton 
energy loss rate (with lepton energy E^) due to ionisation and the losses due 
to Bremsstrahlung, pair-production, hadron production etc. (catastrophic 




(25) 
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losses) is parametrized as 

where P describes the catastrophic loss which dominate over the ionisation 
loss above a certain critical energy ( = a/ (3. This induces a logarithmic 
dependence of the lepton energy loss. 

The calculated ratio TZ varies between ~ 8% to ~ 18% for variation 
from tri-bimaximal mixing scenario and for different values of ^13. Given the 
sensitivity of the ICECUBE detector in terms of detecting GRB neutrino flux 
and considering other uncertainties like that in estimating the flux itself, the 
atmospheric background, low signal yield and the systematic uncertainties of 
the detector, it appears that ICECUBE with its present sensitivity will not 
be able to detect significantly such a small variation due to deviations from 
tri-bimaximal mixing. Hence to detect such small deviation, very precise 
mesurement is called for. This requires more data (more years of run) and 
larger detector size for more statistics. The increase in detector size will not 
widen the deviation of the ratio significantly as the total area factor of the 
detector cancels out in the ratio (Eq. 23) although the total number of both 
muon tracks and total shower yield increase significantly. For the case of 
shower, the whole detector volume is to be considered and from Eq. (22), 
there is indeed an L dependence. This makes the deviation of the ratio wider 
although very marginally as we increase the detector dimension. 

It is difficult to predict the detector dimension and/or the time of expo- 
sure that will be suitable for such a precision measurement discussed above. 
Detailed simulation studies taking into account factors like atmospheric neu- 
trino background, photomultiplier tube efficiency and other possible uncer- 
tainties like the one carried out in Ref. [16] is required for being able to 
comment on the detector parameters for such precise measurements. 

We also want to mention in passing that we have repeated the same 
calculation for single GRBs with fixed red shift {z) values with similar results. 
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Figure Caption 

Fig. 1 Variation of 7^ with 612 and ^23 for (a) ^13 = 0", (b) ^13 = 3", (c) 
^13 = 6" and (d) ^13 = 9°. See text for details. 
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